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ABSTRACT 

It has recently been suggested that the echpsing polar HU Aquarii is host to at least 
two giant planets. We have performed highly detailed dynamical analysis of the orbits 
of those planets and show that the proposed system is highly unstable on timescales 
of < SxlO'^ years. For the coplanar orbits suggested in the discovery letter, we find 
stable orbital solutions for the planetary system only if the outer body moves on an 
orbit that brings it no closer to the host star than ^ 6 AU. The required periastron 
distance for the outer planet lies approximately 5 Hill radii beyond the orbit of the 
inner planet, and well beyond the l-cr error bars placed on the orbit of the outer planet 
in the discovery letter. If the orbits of the proposed planets are significantly inclined 
with respect to one another, the median stability increases slightly, but such systems 
still become destabilised on astronomically minute timescales (typically within a few 
10* years). Only in the highly improbable scenario where the outer planet follows a 
retrograde but coplanar orbit (i.e. inclined by 180° to the orbit of the inner planet) 
is there any significant region of stability within the original l-cr orbital uncertainties. 
Our results suggest that, if there is a second (and potentially, a third planet) in the HU 
Aquarii system, its orbit is dramatically different to that suggested in the discovery 
paper, and that more observations are critically required in order to constrain the 
nature of the suggested orbital bodies. 

Key words: binaries: close, binaries: eclipsing, stars: individual: HU Aqr, planetary 
systems, white dwarfs, methods: N-body simulations 
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d 1 INTRODUCTION 

Precise measurements of timing variations of strictly peri- 
odic events have been successfully used to infer the exis- 
tence bodies orbiting around distant star s. Perhaps the best- 
known examples are the pulsar planets of lWolszczan fc Fraill 
(|l992l ). a system of three extremely low- mass planets or- 
biting the millisecond pulsar PSR1257-I-12. Pulsating sub- 
dwa rf B stars have also b een found to host planets, e.g . V391 
Peg (|Silvotti et al.ll2007^ and HW Vir (|Lee et al.ll2o"ci9i v This 
same technique has recently been applied to eclipsing po- 
lars, using the egress of the small, bright accretion spot as 
a precise "clock." Unseen orbiting bodies can cause small 
shifts in the timing of eclipses (ranging from ~10 sec to 
a few minutes) due to the light-travel time differences im- 
posed by the gravitationa l influence o f the orbiting bodies 
on the system barycentre. iQian et al. discovered the 

first such planet orbiting the eclipsing polar DP Leo by com- 
bining their observed eclipse timings with a long-term data 
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set from ISchwope et all (|2002l ): the combined data set re- 
vealed a sinusoidal variation indicative of a 6.3Mj„p planet 
with a period of 23.8 years. The timing method has shown 
that planets can orbit stars which are wildly different from 
the main-sequence solar-type stars most commonly targeted 
by Doppler and trans it planet search pr ograms. 

In a recent letter, iQian et al.l (|201ll ) announced the dis- 
covery of two (and potentially more) giant planets orbiting 
the eclipsing polar HU Aquarii (hereafter HU Aqr). The 
authors provided fits to the orbits of those planets, plac- 
ing them at orbital radii of 3.6 and 5.4 AU, from the sys- 
tem barycentre, and ascribed minimum masses of 5.9 and 
4.5 Mjup (respectively) to the two bodies. In this Letter, we 
perform a detailed dynamical analysis of the HU Aqr plan- 
etary system in order to assess the stability of the proposed 
planet candidates. 

2 THE HU AQUARII PLANETARY SYSTEM 

In their s t udy of the eclipsing polar system HU Aqr, 
iQian et al.l (|201ll ) consider the temporal variation in the ob- 
served eclipse timings ( "O" ) as compared to predicted tim- 
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Parameter 


HU Aqr (AB)b 


HU Aqr (AB)c 


Eccentricity 


0.0 


0.51±0.15 


Orbital Period (yrs) 


6.54±0.01 


11.96±1.41 


Orbital Radius (AU) 


3.6±0.8 


5.4±0.9 


Minimum Mass (Mj„p) 


5.9±0.6 


4.5±0.5 



Table 1. The orbits of the HU Aqr exoplanets l lQian et al.l201lf) . 



ings ("C") that would be expected from a linear ephemeris. 
By plotting a simple O — C diagram, they show that the 
O — C residuals for HU Aqr contain two cyclical signals 
superposed on a longer-period curvature. Each signal can 
be modelled as a Keplerian orbit to determine the plane- 
tary parameters. We reproduce the parameter estimates of 
iQian et al.l (|201ll ) in Table [T] As for planets detected with 
by radial- velocity method, only the radial component of the 
planet's influence on the host star is detectable. Here, only 
the line-of-sight light-travel time differences are observed, 
so the mass esti mates for the HU A qr planets are given as 
minimum values. iQian et al.l l|201ll ) note that the HU Aqr 
system inclination is 85° , so, if the planets orbit in the same 
plane as the stars, their true masses would only be 0.4% 
larger than the minimum values given in Table [T] 

3 A DYNAMICAL SEARCH FOR STABLE 
ORBITS 

In order to examine the potential dynamical stability of 
the two planets suggested for the HU Aqr system, we per- 
formed a large number of detailed dynamical simulations 
using the Hybrid integrator within the A''-body dynamical 
package Mercury (|Chamberslll999l ). Following the strategy 
employed to analyse t he stability of the HR 8799 system 
(|Marshall. Horner fc Car ter 2010), we held the orbit of the 
inner planet constant (with a = 3.6 AU and e = 0.0), and 
varied the orbital elements of the outer planet across a range 
corresponding to ±3 times the discovery letter's quoted un- 
certainties in the semi-major axis, a, and ecc entricity, e. 
We in itially considered the scenario described in lQian et al.l 
where the planets are considered to be co-planar. In 
other words, we set the orbital inclinations, i, of the two 
planets to be 0° at the start of our integrations. We treat 
the central stars, a 0.88 Mq white dwarf and a 0.2 Mq sec- 
ondary, as a single point mass. Since the stars orbit each 
other with a period of only 2.08 hours (i.e. with a separation 
of 0.004 AU), and the bodies of interest are believed to orbit 
at distances of 3.6 and 5.4 AU, this treatment is dynamically 
ju stified. We give ea ch planet the minimum mass estimated 
in lQian et al.l (|201ll ). We note in passing that, if the masses 
of the prop osed planet s are s ignificantly greater than those 
detailed in iQian et al.l (|201l[ ) , then this could only have a 
deleterious effect on the stability of their proposed orbits. 

Fixing the orbit of the inner planet, we simulated a total 
of 9261 planetary systems. In each simulated system, the in- 
ner planet began on the same orbit, but the orbital elements 
of the outer planet were chosen such that each simulation 
sampled a unique set of possible parameters. We distributed 
the orbital elements of the outer planet such that we tested 
21 values of semi-major axis, sp read evenly across =fc3-a from 
the value (a = 5.4 AU) given in lQian et al.l l|201ll ). For each 
value of semi-major axis, we tested 21 values of orbital ec- 



centricity, again spread evenly across ±3-(t from the value 
(e = 0.51), given in that work. Finally, at each of these 
441 (a, e) locations we carried out 21 tests, with the initial 
location of the planet distributed across a range of ±3-cr 
from the nominal m ean anomal y of th e outer planet (calcu- 
lated from Fig. 2 of iQian et all (|201ll ')'). These 9261 unique 
models, based on the HU Aqr system parameters, were in- 
tegrated using the Hybrid integrator within Mercury for a 
period of 100 Myr, following the evolution and final fates 
of the two postulated planets in the system. A planet was 
deemed ejected from the system upon reaching a distance 
of 1000 AU from the barycentre, and all mutual collisions 
were recorded. This yielded a lifetime from each individual 
integration in the range - 100 Myr, defined as the time un- 
til one or other of the planets was removed from the system 
through either collision or ejection. 

To explore the 1-a parameter range in greater detail, 
we launched a second suite of integrations, again using 9261 
test systems. The set-up was performed exactly as described 
above, except that the orbital parameters of the outer planet 
were varied within a 1-a range, rather than the 3-(t distri- 
bution carried out previously. 

These two suites of integrations yielded 17493 distinct 
tests of the stability of the HU Aqr system, with 9261 of 
these perf ormed in the centr al ±l-cr of the element space de- 
scribed in iQian et all (|201ll ). and the other 8232 distributed 
in the range 1-cr to 3-cr. The results of these integrations are 
shown in Fig[T] 

Once these simulations had been completed, we carried 
out equivalent suites for scenarios where the outer planet 
was moving on an orbit inclined to that of the inner planet. 
Following the procedure detailed above, we tested systems 
in which the outer planet's orbit was inclined by 5°, 15° 
and 45° with respect to the inner's orbit, in order to exam- 
ine the influence of mutual orbital inclinations on the sta- 
bility of the system. We then considered further scenarios 
in which the outer planet was moving in a retrograde sense, 
with respect to the inner body, with inclinations of 135° and 
180° . For simplicity, we kept the masses of the two planets 
co nstant thro u gh th ese runs at the lowest values suggested 
m iQian et all (|201ll ). Although it is true that significantly 
inclined orbits for the planets would result in larger real 
masses for them, we note that a significant mutual inclina- 
tion between the planets does not necessarily mean that it 
is the outermost planet that is inclined to the line of sight, 
while the innermost is in that plane. Rather than attempt 
to shift the planetary mass by the free parameter of poten- 
tial inclination, we instead took the lowest masses possible. 
We remind the reader that this essentially means that we 
have allowed the planetary system the best possible chance 
of being stable - as the mass of the planets increases, so does 
their gravitational reach, increasing the strength of any mu- 
tual interactions. 

4 RESULTS 

The results of our i =0° (i.e. coplanar) dynamical integra- 
tions are shown in Fig [T] Each colour box in that figure 
shows the median lifetime obtained from 21 independent 
integrations performed with the outer planet placed on an 
orbit with that particular combination of a and e. The lo- 
cation of the inner planet is marked with a hollow circle, 
whilst the location of the nominal orbit for the outer planet 
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given in lQian et al.l l|201lh is shown by a hollow square, with 
the 1-cr error bars given in that work denoted by the solid 
lines stretching from that square. The vertical dot-dash lines 
show the locations of the strongest mean-motion resonances 
(hereafter MMR) with the orbit of the inner-most planet. 
The curved dotted lines that meet at the location of the in- 
ner planet join all orbits whose periastron (outward curving 
line) or apastron (inward curving line) lie at a distance equal 
to the orbital radius of the inner planet. Every point in the 
region above these lines denotes orbits for the outer planet 
that cross that of the inner. The two dotted lines labelled 
"ARhui and SRhui connect orbits which pass periastron at a 
distance equal to the orbital radius of the inner planet plus 
3 and 5 times that planet's Hill radius (Rh), respectively, 
where Rh is defined as 



Rh = cip 



1/3 



(1) 



The Hill radius, Rh, is commonly used in studies of 
orbital dyna mics as a proxy for the dy namical "reach" of a 
given body (|Horner et all 12001 l2004al lb[). Close encounters 
between two massive bodies are typically defined as those 
that occur at a distance closer than 3 Rh, although some 
particularly conservative studies consider 5 Rh a sufficiently 
close approach to be labeled as such. These lines therefore 
show the limits at which the outer planet approaches the 
inner within the prescribed number of Rh - those orbits 
outside the lines are too widely separated at the start of 
the integrations to undergo close encounters, while all those 
within the region bounded by the periastron and apastron 
lines for the orbital radius of the inner planet are orbits 
which cross that of that planet. 

It is immediately apparent in Fi g [1] that the great ma- 
jority of possible orbits suggested in Qia n et al.l ( 201 ll ) are 



extremely unstable, with just a small region at low eccentric- 
ities and high semi-major axes (i.e. below the 3 Rh line in the 
bottom right-hand corner) showing any significant long term 
stability. This is not a surprising result - the two planets in 
question have particularly large masses, and therefore have 
very large dynamical reaches - and so must be widely sepa- 
rated in order that they do not strongly perturb one another. 
Indeed, it is clear from that figure that the regions of stabil- 
ity and instability for the HU Aqr system are a strong func- 
tion of the peri- and apastron distances of the initial orbits 
of the outer planet. The only orbits that display relatively 
strong stability all have periastron at distances greater than 
5 Rh beyond the orbit of the inner planet. Any closer, and 
the encounters between the planets are strongly disruptive. 
Orbits of the outer planet which approach the inner planet 
to a distance between 3 and 5 Rh are clearly more stable 
than those which come closer to that planet, but still display 
significant instability on astronomically short timescales. 

However, what of mutual inclinations between the two 
planets? Could the system be stabilised by the planets mov- 
ing on orbits that are significantly inclined with respect to 
one another? This question can be answered by examina- 
tion of Fig[2l which presents the results for all six scenarios 
considered in this work. 

As the orbital inclination of the outer planet is increased 
to 5° and then 15°, the overall stability of the system ap- 
pears to increase somewhat, with the yellow and orange 
colours that denote moderate stability spreading across the 
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Figure 1. Lifeti me stability plot of the HU Aqr system using 
the parameters of lQian et al.l {ioi^) with the planets on coplanar 
orbits. The simulations detailed cover the 3-(t parameter space 
in a grid of 21x21 points. The central 1-cr parameter space is 
covered by an additional and denser set of 21x21 points. Each 
grid point is the median lifetime of 21 simulated HU Aqr systems 
with the outer planet's initial a and e. The location of the inner 
planet is denoted by the point of the solid triangle on the x-axis, 
while the nominal orbit of the outer planet is marked by the 
filled square. Solid lines show the ext ent of the l-cr errors in a 
and e suggested by lOian et al.l l|201ll ') . The vertical dot-dashed 
lines show the location of the strongest MMRs in relation to the 
orbit of the inner planet. The two dotted lines radiating from 
the location of the inner planet connect all orbits of that have 
either their periastron (outward curving line) or apastron (inward 
curving line) at the location of the planet. As such, all orbits in the 
region bounded by these lines cross the orbit of the inner planet. 
The dotted lines labelled 3 Rh and 5 Rh connect all orbits that 
pass periastron at a distance of three or five Hill radii beyond the 
orbit of the inner planet. It can readily be seen that the planetary 
system shows extreme instability, aside from in a small region of 
a-e space to the lower right hand side of the figure, corresponding 
to orbits of the outer planet that remain at a barycentric distance 
beyond ~6 AU. 



entire plot. However, the effect of increased inclination on 
the region of greatest stability, those orbits with median life- 
times greater than one million years, is negligible. Indeed, 
by the time the inclination is increased to 15°, that region 
appears to shrink somewhat, with particular destabilisation 
occurring just inwards of the location of the 5:2 MMR with 
the orbit of the inner planet. Once the orbital inclination of 
the outer planet is increased to 45° (top right hand panel), 
this effect becomes far more pronounced, with only small re- 
gions of stability remaining in low-eccentricity orbits around 
the 2:1 MMR, around the 5:2 MMR, and a larger region be- 
yond the location of the 3:1 MMR. Indeed, far from enha 
ncing the stability of the planetary system, as might be ex- 
pected, the 3:1 MMR acts to destabilise the orbits of the 
planets in this more highly inclined scenario. 

The central panel on the right hand side of Fig [5] shows 
the scenario for which the orbit of the outer planet is in- 
clined by 135° to that of the inner. That scenario displays 
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Figure 2. Lifetime a-e stability plots of the HU Aqr system using the parameters of iQian et al ] l|201ll l. The top-left panel reproduces 
the results shown in Fig[T] detailing the scenario in which both planets are considered to be moving on co-planar, prograde orbits. The 
panel to the middle-left shows the results for the scenario where the outer planet has an orbit inclined by 5° to that of the inner planet, 
with the lower-left, top-right, middle-right and lower-right panels showing the results for initial orbital inclinations of 15°, 45°, 135° and 
180°, respectively. The various features plotted on each panel are the same as those shown in Fig[T] 



a remarkable lack of stability across the entire 3-a region of 
plausible orbits for the outer planet. 

The only scenario in which orbits within the 1-a un- 
certainty range for the outer planet display strong stability 
is shown in the lower-right hand panel of Fig. [5] (i =180°), 
when the orbit is both retrograde and coplanar. In that sce- 
nario, almost all orbits for the outer planet not crossing 
that of the inner are dynamically stable on long timescales. 
In that extreme scenario, even some configurations in which 
the orbits of the two planets cross show significant stabil- 
ity. Those scenarios are the only ones of our entire suite of 
almost 105,000 test integrations in which orbits of the HU 
A qr planets si t uated within the 1-a error bounds stipulated 
m iQian et all (|201lD can display any long-term stability. 

5 CONCLUSIONS AND DISCUSSION 

We have investigated the stability of the recently discovered 
planetary system around HU Aqr across the ±3-(t uncer- 
tainty ranges for the orbit of the outer planet. In the sim- 
plest, coplanar, case, we find that the system is unstable on 
timescales < 5 x 10"^ years, except for scenarios in which the 
mutual separation between the two planets was greater than 



5 Rh at the outer planet's periastron. When the orbit of the 
outer planet is set such that it reaches periastron between 
3 Rh and 5 Rh beyond the orbit of the inner planet, the or- 
bital stability is somewhat greater than scenarios where the 
minimum separation between the two planets is smaller then 
3 Rh, but the planetary system still falls apart on timescales 
too short to inspire any confidence that the tested orbits rep- 
resent the true state of the HU Aqr system. We see a steady 
increase in the median lifetime of the system across much of 
the tested phase-space as the mutual inclination of the orbits 
climbs from the coplanar case, through 5° to 15° to 45°, as 
expected from the resulting reduction in the amount of time 
the two planets spend in close proximity. However, even in 
the most extreme prograde case investigated (i = 45°), the 
lifetime of systems within the 1-a parameter space did not 
exceed 10^ years. This is much shorter than the expected 
system age, suggesting that we have either caught the sys- 
tem during a period of significant dynamical instability in 
which the planets are undergoing a rearr angement /e jection , 
or that the orbital solution described in iQian etal (|201ll ) 
does not represent the true state of the system. Interest- 
ingly, although an increase in the mutual inclination of the 
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planetary orbits causes those orbits in the unstable region 
to become slightly more stable, it also results in a reduction 
in the stable region at distances greater than 5 Rh from the 
inner planet. 

If the orbit of the outer planet is retrograde compared 
to that of the inner planet, but the orbits remain coplanar 
(i.e. the mutual inclination of the two orbits is 180°, then 
we find that the system could be stable across a wide range 
of parameter space, including some scenarios within the 1- 
(j errors quoted for the orbit of the outer planet. However, 
it seems difficult to comprehend how the proposed planets 
could have evolved into such orbits. By contrast, if the orbits 
of the two planets have mutual inclinations of 135°, then we 
find that no region of the tested a-e phase space is stable on 
timescales greater than ~ 10* years. 

The results of our dynamical simulations raise a num- 
ber of interesting possibilities . Assu ming that the detection 
of the planets bv lQian et al.l l|201ll ) is robust, and that the 
orbits of the planets at the current epoch are exactly as de- 
scribed in that work, the planetary system must be going 
through a dramatic period of dynamical instability, which 
will in short order result in the loss of one (or both) of the 
planets therein. However, the incredibly short lifetimes we 
find for the proposed planets suggest that, statistically, this 
is unlikely (the odds of observing a planetary system during 
the last few thousand years of a multi-billion-year lifetime 
seem remarkably small). On the other hand, if we assume 
that the detection of the planets is robust, but the orbital 
parameters given are not a good measure of the true state of 
the system, we suggest that it is most likely that the outer 
planet is moving on a low-eccentricity orbit far from the cen- 
tral bodies. Such an orbit would allow that planet to remain 
sufficiently far from the inner planet to be dynamically sta- 
ble on multi-million year timescales, and therefore seems a 
much more reasonable solution for t he dynamics of th e sys- 
tem. If the two planets detected bv lQian et all (120111 ) have 
significant mutual orbital inclinations, then we find that a 
slightly wider range of orbits are possible for the outer planet 
to display moderate stability, but that at the same time the 
region of greatest stability for that planet's orbit actually 
reduces in size. In any case, in such scenarios the stable re- 
gions remain relatively restrictive both in terms of orbital 
eccentricity and semi-major axis - keeping the planet well 
beyond the inner's sphere of influence. 

The only way in which the orbits of the two planets 
can be induced to show significant stability is to consider a 
scenario in which they are coplanar, but with the orbit of 
the outer planet being retrograde with respect to that of the 
inner. In such a scenario, a wide area of the studied a-e phase 
space becomes dynamically stable, inclu ding some solution s 
within the 1-a uncertainties detailed bv lQian et al.l (|201ll) . 
Whilst this might appear promising, we note that it is hard 
to envision a way in which the planets could evolve into such 
an unusual configuration without significantly destabilising 
one another's orbits. 

When one examines the residuals from the 0-C dia- 
gram in the lQian et al.l (|201ll ) work, a first glance suggests 
that a low-eccentricity orbit for the outer planet is incom- 
patible with the observed data. However, it is p ossible that 
the re moval of the long-term quadratic trend bv lQian et al.l 
l|201ll ) could have resulted in an artificially enhanced ec- 
centricity for the outer planet, reducing its stability. On the 



other hand, detailed simulati ons of Doppler velocity data by 
lAnglada-Escude et al.l l|2010l ) demonstrated that two plan- 
ets in circular orbits can mimic the signal of a single planet 
on an eccentric orbit, so long as those planets move on mutu- 
all y resonant orbits. T his possibility was recently explored 
by iTinnev et ai] l|201ll ) for the case of HD 38283b, an ec- 
centric planet in a one- year orbit. If, rather than a highly 
eccentric second planet, we have a scenario where the sys- 
tem contains at least three massive planets, each moving on 
low-eccentricity orbits, it might be poss ible to expl a in the 
suggested shape of the O-C diagram in lOian et al.l (1201 ll ) 
whilst placing the planets on or bits that are dynamically 
stable. Indeed, lOian et al.l (|201l[ ) suggest that there might 
be a third massive body at a large barycentric distance, 
based on the archaic Titius-Bode law. While the use of that 
law is not generally encouraged as a predictive tool in exo- 
planetary science, the presence of a more distant third body 
would allow orbital fits with the second planet moving on 
a low-eccentricity orbit. Although we have not dynamically 
simulated such a speculative scenario, it seems reasonable 
to assume that if a third planet were located at an orbital 
radius at least 5 Rh beyond that of the second planet, the 
system could display long term dynamical stability without 
violating the observed variations in the egress of the accre- 
tion hotspot from eclipse by the stellar secondary body. 

A more detailed statistical analysis of this highly fasci- 
nating exoplanetary system is clearly necessary in order to 
disentangle the true nature of the proposed planetary sys- 
tem. Such work will doubtless throw fresh light on one of 
the most peculiar planetary systems detected to date. 

ACKNOWLEDGMENTS 

JH gratefully acknowledges the financial support of the Aus- 
tralian government through ARC Grant DP0774000. RW is 
supported by a UNSW Vice-Chancellor's Fellowship. JPM 
is partly supported by Spanish grant AYA 2008/01727, and 
gratefully acknowledges Maria Cunningham for funding his 
collaborative visit to UNSW. We also wish to thank the 
anonymous referee for their swift and very helpful feedback. 

REFERENCES 

Anglada-Escude, C, Lopez-Morales, M., & Chambers, 

J. E. 2010, ApJ, 709, 168 
Chambers, J. E., 1999, MNRAS, 304, 793 
Horner, J., Evans, N. W., Bailey, M. E., & Asher, D. J. 

2003, MNRAS, 343, 1057 
Horner, J., Evans, N. W., & Bailey, M. E. 2004, MNRAS, 

354, 798 

Horner, J., Evans, N. W., & Bailey, M. E. 2004, MNRAS, 

355, 321 

Lee, J. W. et al., 2009, AJ, 137, 3181 

MarshaU, J. P., Horner, J. & Carter, A., 2010, IJA, 9, 259 
Qian, S.-B., Liao, W.-P., Zhu, L.-Y., & Dai, Z.-B. 2010, 

ApJL, 708, L66 
Qian, S.-B. et al., 2011, MNRAS, doi:10.1111/j.l745- 

3933.2011.01045.x 
Schwope, A. D., Hambaryan, V., Schwarz, R., Kanbach, 

C, Gansicke, B. T. 2002, A&A, 392, 541 
Silvotti, R., et al. 2007, Nature, 449, 189 
Tinney, C. G. et al., 2011, ApJ, 732, 31 
Wolszczan, A., & FraU, D. A. 1992, Nature, 355, 145 



